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Abstract. Three closely related stumbling blocks of solar mean field dynamo theory are discussed: how dominant are the 
small scale fields, how is the alpha effect quenched, and whether magnetic and current helicity fluxes alleviate the quenching? 
It is shown that even at the largest currently available resolution there is no clear evidence of power law scaling of the magnetic 
and kinetic energy spectra in turbulence. However, using subgrid scale modeling, some indications of asymptotic equipartition 
can be found. The frequently used first order smoothing approach to calculate the alpha effect and other transport coefficients 
is contrasted with the superior minimal tau approximation. The possibility of catastrophic alpha quenching is discussed as a 
result of magnetic helicity conservation. Magnetic and current helicity fluxes are shown to alleviate catastrophic quenching 
in the presence of shear. Evidence for strong large scale dynamo action, even in the absence of helicity in the forcing, is 
presented. 
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1. Introduction 

Over the past 30 years, the standard approach to understand- 
ing the origin of the solar cycle has been mean field dynamo 
theory. This approach can be justified simply by the fact that 
the sun does have a finite azimuthally averaged mean field, 
B(r,9,t), where r and 9 are radius and colatitude. Obser- 
vationally, we really only know with some certainty its ra- 
dial component at the surface, B r (R,8,t), where R is the 
radius of the sun. There is also some indirect evidence for 
the toroidal field, -B0(r spo t, 8, t), where r spo t is the not well 
known radius where sunspots are anchored. Both components 
give a very clear indication of the spatio-temporal coherence 
of the mean field, with a 22 year cycle and latitudinal migra- 
tion. 

Mean field theory has certainly been successful in show- 
ing that a solar-like mean field can be produced if there is an 
a effect, i.e. if the mean electromotive force has a component 
along the mean field (Weiss 2005). A number of complica- 
tions have arisen in the mean time. 
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(i) The standard theory for calculating the value of a and 
other relevant transport coefficients (such as turbulent 
magnetic diffusivity) relies on the first order smooth- 
ing (FOSA) or second order correlation approximation 
(SOCA). This approach is valid if either the magnetic 
Reynolds number is small (poor microscopic conductiv- 
ity of the gas) or if the so-called Strouhal number (Kapyla 
et al. 2005) is small compared to unity. The latter means 
that the correlation time is supposed to be much less than 
the turnover time of the turbulence, which is not usually 
the case. The former assumption is equally inappropriate. 

(ii) Shortly after mean field theory became popular there have 
been recurrent concerns about its applicability when the 
field strength is comparable to the equipartition value of 
the turbulence, i.e. if the magnetic energy is compara- 
ble to the kinetic energy of the turbulence. Such con- 
cerns where first addressed by Piddington (1970, 1972), 
but more recently by Vainshtein & Cattaneo (1992) and 
Kulsrud & Anderson (1993). 

(iii) In more recent years the problem of the resistively slow 
evolution of magnetic helicity has been discussed in con- 
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nection with a correspondingly slow saturation if the dy- 
namo (Brandenburg 2001, Mininni et al. 2003). This re- 
sistively slow time scale in the problem can also af- 
fect the cycle period in an ail dynamo (Brandenburg et 
al. 2001, 2002). The presence of boundaries may help 
(Blackman & Field 2000a,b, Kleeorin et al. 2000, 2002, 
2003), but it may also make matters worse (Brandenburg 
& Dobler 2001). The question is therefore whether suit- 
ably arranged shear is needed to transport magnetic he- 
licity out of the domain (Vishniac & Cho 2001, Subrama- 
nian & Brandenburg 2004). The indications are that with 
shear and open boundary conditions the otherwise cata- 
strophic a quenching can be alleviated (Brandenburg & 
Sandin 2004). 

In the following we discuss the current status on all three 
issues. Indeed, there has been a lot of progress, and from an 
optimistic viewpoint one might almost think that these issues 
have now been solved. But, of course, as always in science, 
new problems and unexpected issues emerge all the time. 
Also, having proposed one solution to a problem does not 
exclude alternative solutions, and so only with time will we 
be able to look back and say how it really was. 



2. Problem I: dominance of small scale fields? 

The first issue of the relative importance of small scale versus 
large scale fields can be discussed in terms of the magnetic 
energy spectrum. We begin with a historical perspective. 

2.1. Turbulent diffusion and turbulent cascade 

Turbulent diffusion relies on the ability of the turbulence to 
transport energy from large scales to small scales. This is an 
important consideration that was discussed by Stix (1974) in 
the context of early criticism raised against dynamo theory. 
The situation is illustrated in Fig.Q 

The assumption of a magnetic energy spectrum being par- 
allel or even coincident with that of kinetic energy, and this 
being equivalent to the cascade in ordinary (nonmagnetic) 
turbulence is natural, but not trivial and not fully confirmed 
by simulations, as will be discussed in the next section. 

2.2. Numerical indications for a magnetic cascade 

Over the past decades some steady progress has been brought 
about; see Fig. 13 where we show kinetic and magnetic en- 
ergy spectra from Meneguzzi et al. (1981), Kida et al. (1991), 
Maron & Cowley (2001), representing the improvement of 
the state of the art simulations over the past three decades. 
The spectra from all three papers show that at large scales 
the spectral magnetic energy is below the spectral kinetic en- 
ergy. In the spectra of Meneguzzi et al. (1981) and Maron & 
Cowley (2001) the spectral magnetic energy peaks at about 
k m 5. However, there has not really been any evidence for 
power law behavior. One exception is the spectrum obtained 
by Kida et al. (1991) who find a fc° spectrum, but this has 




Fig. 1. A sketch from the original review by Stix (1974) 
where he illustrates the role of a turbulent cascade in con- 
verting large scale poloidal and toroidal flux into small scale 
fields where it can be converted into heat by microscopic re- 
sistivity. 



not been found in subsequent simulations by Maron & Cow- 
ley (2001) or Maron et al. (2004), for example. The numeri- 
cal resolution used in the three cases is 64 3 , 128 3 , and 256 3 
meshpoints, respectively. Obviously, much higher resolution 
is needed to begin to address the possibility of powerlaw be- 
havior. This was the main reason for Haugen et al. (2003) to 
push the resolution to 1024 3 meshpoints; see Fig. [5] Based on 
these results, one can begin to see the development of what 
looks like an inertial range with a tentative /c~ 3//2 scaling. 
If this is confirmed, it would rule out earlier claims that the 
magnetic energy spectrum peaks at the resistive scale (Maron 
& Blackman 2002, Schekochihin et al. 2002). However, the 
results by Haugen et al. (2003) are still open to alternative in- 
terpretations. Schekochihin et al. (2004a) have argued that the 
magnetic spectrum is still curved, and that there is therefore 
actually no evidence for powerlaw behavior. So, the asymp- 
totic spectral behavior of hydromagnetic turbulence is still 
very much an open question. 

A puzzling aspect of all the magnetic energy spectra is 
the excess of spectral magnetic energy over spectral kinetic 
energy. In order to get some insight into the possible asymp- 
totic behavior, Haugen & Brandenburg (2004) have recently 
considered simulations with subgrid scale modeling using ei- 
ther hyperviscosity or Smagorinsky viscosity together with 
magnetic hyper-resistivity. 
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Fig. 2. Steady progress in solving the small scale dynamo problem. Kinetic and magnetic energy spectra from the papers by 
Meneguzzi et al. (1981), Kida et al. (1991), and Maron & Cowley (2001). The numerical resolution used in the three cases 
is 64 3 , 128 3 , and 256 3 meshpoints, respectively. All spectra agree in that at large scales (small wavenumbers k) the spectral 
magnetic energy is below the spectral kinetic energy. In the spectra of Meneguzzi et al. (1981) and Maron & Cowley (2001) 
the spectral magnetic energy peaks at about k = 5. 
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Fig. 3. Magnetic, kinetic and total energy spectra. 1024 3 
meshpoints. The Reynolds number is w rms /(yfcf) ~ 960. 
[Adapted from Haugen et al. (2003).] 
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Fig. 4. Comparison of energy spectra of the 4096 3 mesh 
points run of Kaneda et al. (solid line) and 512 3 mesh points 
runs with hyperviscosity (dash-dotted line) and Smagorinsky 
viscosity (dashed line). 



Consider for comparison first the purely hydrodynamic 
case (Fig. |4}, where we can use the high resolution simula- 
tions of Kaneda et al. (2003) on the Earth Simulator as 
benchmark. This simulation corresponds to a resolution of 
4096 collocation points. Note that there is not even a clear 
confirmation of the famous Kolmogorov fc -5 / 3 spectrum, but 
there is a 1 correction in what looks like the inertial 
range. However, simulations with both hyperviscosity (Hau- 
gen & Brandenburg 2004) and Smagorinsky subgrid scale 
modeling (Haugen & Brandenburg, in preparation) confirm 
this correction to the inertial range; see Fig.|4] 

The uprise in the compensated power spectra just before 
the dissipative subrange is due to the bottleneck effect in tur- 
bulence (Falkovich 1994). This effect is much weaker in wind 
tunnel turbulence (She & Jackson 1993), but this is because 
these are one-dimensional spectra, E^q (k), which are related 



to the fully three-dimensional spectra via a simple integral 
transformation 

EMk)=f^Pdk>. (1) 

Of course, for perfect power law spectra the two are the same, 
but they can be quite different when there are departures from 
power law behavior, such as due to the bottleneck effect itself 
and due to the dissipative subrange (Dobler et al. 2003). 

In Fig.|5]we show the results for the hydromagnetic case 
where we have used hyperresistivity and either hyperviscos- 
ity or Smagorinsky subgrid scale modeling for the velocity 
field. There are two important things to notice. First, the com- 
pensated magnetic spectrum seems flat, but the kinetic energy 
spectrum seems to rise, possibly approaching the magnetic 
energy spectrum. A sketch of what we believe the asymptotic 
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Fig. 5. Magnetic and kinetic energy spectra for runs with 
512 3 mesh points and hyperviscosity and hyper-resistivity 
(solid line) and Smagorinsky viscosity and hyper-resistivity 
(dashed line). Note the mutual approach of kinetic and mag- 
netic energy spectra before entering the dissipative subrange. 
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Fig. 6. Sketch of the anticipated kinetic and magnetic energy 
spectra in the large Reynolds number limit for Pr m = l.Note 
the slight super-equipartition just to the right of the peak of 
Em (k) and the asymptotic equipartition for large wavenum- 
bers. 

magnetic and kinetic energy spectra could look like is shown 
in Fig.|6] 

In summary, there seems now some evidence suggesting 
that there is indeed a magnetic energy cascade. We should 
emphasize, however, that most of the simulations to date have 
either unit magnetic Prandtl number, or at least magnetic 
Prandtl numbers that are not very different from unity. 

The question of the magnetic Prandtl number dependence 
is potentially quite important for stars where this number is 
very small (magnetic diffusivity much larger than the kine- 
matic viscosity). Using a modified Kazantsev model, Ro- 
gachevskii & Kleeorin (1997) pointed out that when the ve- 
locity field displays Kolmogorov scaling near the resistive 
cutoff wave number, the dynamo becomes much harder to 
excite than for a velocity field that shows significant power 
only at large scales. They found that the critical magnetic 
Reynolds number for dynamo action can exceed a value of 
400, which is about 10 times larger than the critical value 
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Fig. 7. Dependence of the critical magnetic Reynolds num- 
bers as a function of the magnetic Prandtl number using for 
the direct simulations (Haugen et al. 2004), indicated by di- 
amonds, compared with results using hyperviscosity (aster- 
isks) or Smagorinsky viscosity (triangles). 

of 35 for unit magnetic Prandtl number. A similar result has 
been found more recently by Schekochihin et al. (2004b) and 
Boldyrev & Cattaneo (2004). Direct simulations suggest that 
the critical magnetic Reynolds number increases with de- 
creasing magnetic Prandtl number, at least like i? mjC rit ~ 

— 1/2 

35 Pr M (Haugen et al. 2004); this result is not really con- 
firmed by simulations with Smagorinsky viscosity and hy- 
perviscosity which seem to suggest a sharper dependence 
(Fig. 0; see Schekochihin et al. (2004c). One might have 
thought that this could be connected with the artificially en- 
hanced bottleneck effect in the hyperviscous simulations, but 
the idea is not supported by the direct inspection of the kinetic 
energy spectra that do not show a bottleneck effect. 

In simulations with helicity, that will be discussed in a 
later section, the value of the magnetic Prandtl number does 
not seem to affect the onset of dynamo action. Indeed, simu- 
lations with a magnetic Prandtl number of 0. 1 (Brandenburg 
2001) showed that the critical value of the magnetic Reynolds 
number was the same as for unit magnetic Prandtl number. 

3. Problem II: is first order smoothing OK? 

Let us now turn to the question of how to generate large scale 
fields. This question can be addressed in terms of mean field 
dynamo theory where the a effect (and perhaps other effects) 
play a central role. The value of a is traditionally calculated 
using the first order smoothing approximation, where one ne- 
glects triple correlations. This approach is in principle not 
applicable when the magnetic Reynolds number is large, but 
it seems to work anyway. Why this is so, and how one can do 
better, can perhaps best be understood in the simpler case of 
passive scalar diffusion. 

3.1. Passive scalar diffusion 

In this section we describe both the first order smoothing ap- 
proximation and the minimal tau approximation. We follow 
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3. 1 Passive scalar diffusion 



here the presentation of Blackman & Field (2003), who stud- 
ied the passive scalar case as a simpler test case of the more 
interesting magnetic case which they studied earlier (Black- 
man & Field 2002). The minimal tau approximation was first 
used by Vainshtein & Kitchatinov (1986) and Kleeorin et al. 
(1990). 

The evolution equation for the passive scalar concentra- 
tion C is, in the absence of microscopic diffusion, 

DC - (2, 



Df 



= 



where D/Di = d/dt + u ■ V is the advective derivative. 
We assume, for simplicity, that the flow is incompressible, 
i.e. V ■ u — 0, and that C shows some large scale variation 
so that a meaningful average can be defined (denoted by an 
overbar). Thus, after averaging we have 



dC 
dt 



d 

oxj 



(3) 



One is obviously interested in a closed equation for the mean 
flux of passive scalar concentration, 



T — uc. 



(4) 



in terms of the mean concentration, i.e. T = J~(C). For the 
reader who is more familiar with the magnetic case we note 
that one should think of T being similar to the electromotive 
force £ = u x b where, in turn, one is interested in a closed 
equation of the form £ = £(B). 

Subtracting Eq. (|3j from Eq. (0 we obtain the equation 
for the fluctuation c = C — C, 

dc — 

— = -djiujC) - d 3 (u jC ) - djUje. (5) 

This is where we come to a turning point. In the first order 
smoothing approximation (FOSA) one calculates 



(6) 



F{t) = u{t) J c(t') dt' (FOSA), 

where c = dc/dt, and we have omitted the common x de- 
pendence on all quantities. In the minimal tau approximation 
(MTA) one calculates instead 



dT — — 

-— =uc + iic (MTA). 

ot 



(7) 



A principle difference between the two approaches is that the 
momentum equation is naturally incorporated under MTA, 
i.e. 



ii = —UjdjUi 



(8) 



where p is the pressure. A more interesting situation would 
arise if we allowed here for rotation and included the Coriolis 
force, but we omit this here. Thus, we have 



Qj: _ 

— — = —UiUj dC —UiUjdjC — cujdjUi — cdip . 
ot ■> , ' 

triple correlations =T 



(9) 



This equation shows the important point that, at least in the 
steady state, the triple correlations are never negligible. In- 
stead, they are actually comparable to the quadratic correla- 
tion term on the right hand side. (We remark that the first two 



terms in T cancel if the volume averages are used to inte- 
grate by parts, but the third term still remains.) The closure 
hypothesis used in MTA states that 



T = (MTA closure hypothesis). (10) 

r 

Inserting this expression and moving this T to the left hand 
side, we have 



dt 



(11) 



The extra time derivative on the right hand side corresponds 
to the Faraday displacement current in electrodynamics. In 
the present case, this term can be neglected if the diffusion 
speed is less than the turbulent rms velocity. In some sense 
this is of course never the case, because with ordinary Fick- 
ian diffusion (time derivative neglected) the diffusion process 
is described by an elliptic equation which has infinite signal 
propagation speed, and is hence violating causality. 

In the following section we discuss recent work by Bran- 
denburg et al. (2004) who showed that the presence of this 
displacement term in a non-Fickian version of the diffusion 
equation is indeed justified. Before turning to this aspect, let 
us contrast the MTA closure with the FOSA closure assump- 
tion. Here one is instead dealing with an integral equation, 



T,. = 



Ui(t)uj (f) 8C(t') dt' + triple correlations, (12) 



where the triple correlations are neglected [unless higher or- 
der terms are included; see Nicklaus & Stix (1988) or Car- 
valho (1992)]. 

In this particular case the two approaches become equiv- 
alent if one assumes that the two-times correlation function 
Ui(t)uj(t') is proportional to muj exp[— (t — f)/r]. This 
can be shown by differentiating Eq. ill 2b . Thus, the causality 
problem in the Fickian diffusion approximation stems really 
only from the commonly used approximation that the two- 
times correlation function can be approximated by a delta 
function. If the extra time derivative is not neglected, the dif- 
fusion equation becomes a damped wave equation, 



d 2 C 1 dC 



dt 2 



— I,,* \7 z r i 

Q-f- ~ 3 "rms v ^ ' 



(13) 



where the wave speed is u rms /\/3. Note also that, after mul- 
tiplication with r, the coefficient on the right hand side be- 
comes ^Tu 2 ms = n t , and the second time derivative on the 
left hand side becomes unimportant in the limit r — > 0, or 
when the physical time scales are long compared with r. In 
that case we have simply 



dC_ 

dt 



K t V 2 C. 



(14) 



In the following we discuss a case where Eq. dl4> is clearly 
insufficient and where the full time dependence has to be re- 
tained. 



3.3 Significance for the magnetic case 
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Fig. 8. Passive scalar amplitude, (C ) 1 ^ 2 , versus time (nor- 
malized by u rms fcf) for two different values of kf/ki. The 
simulations have 256 3 meshpoints. The results are compared 
with solutions to the non-Fickian diffusion model. [Adapted 
from Brandenburg et al. (2004).] 
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Fig. 9. Strouhal number as a function of kf/k\ for differ- 
ent values of ReLs, i e- the large scale Reynolds number. The 
resolution varies between 64 3 meshpoints (ReLs = 100) and 
512 3 meshpoints (Re LS = 1000). 

3.2. Turbulent displacement flux and value of r 

A particularly obvious way of demonstrating the presence of 
the second time derivative is by considering a numerical ex- 
periment where C = initially. Equation MAI would predict 
that then C — at all times. But, according to the alternative 
formulation (1131 . this need not be true if initially dC/dt ^ 0. 
In practice, this can be achieved by arranging the initial fluc- 
tuations of c such that they correlate with u z . Of course, such 
highly correlated arrangement will soon disappear and hence 
there will be no turbulent flux in the long time limit. Never- 
theless, at early times, (C ) 1 / 2 (an easily accessible measure 
of the passive scalar amplitude) rises from zero to a finite 
value; see Fig. [8] 

Closer inspection of Fig. [8] reveals that, when the 
wavenumber of the forcing is sufficiently small (i.e. the 
size of the turbulent eddies is comparable to the box size), 
(C ) x l 2 approaches zero in an oscillatory fashion. This re- 



markable result can only be explained by the presence of the 
second time derivative term giving rise to wave-like behavior. 
This shows that the presence of the new term is actually justi- 
fied. Comparison with model calculations shows that the non- 
dimensional measure of r, St = ru lms kf, must be around 

3. (In mean-field theory this number is usually called the 
Strouhal number.) This rules out the validity of the quasilin- 
ear (first order smoothing) approximation which would only 
be valid for St — > 0. 

Next, we consider an experiment to establish directly the 
value of St. We do this by imposing a passive scalar gradi- 
ent, which leads to a steady state, and measuring the result- 
ing turbulent passive scalar flux. By comparing double and 
triple moments we can measure St quite accurately without 
invoking a fitting procedure as in the previous experiment. 
The result is shown in Fig. |9] and it confirms that St w 3 in 
the limit of small forcing wavenumber, k{ . The details can be 
found in Brandenburg et al. (2004). 

3.3. Significance for the magnetic case 

In the hydromagnetic case one has £ = uxb + uxb. Here 
the u term is particularly important and usually not obtained 
with FOSA. Focusing on the term Bibi.i on the right hand 
side of the Uj equation, we get 

p (ii x b)i = eijkb k b jt p B p + ... l = \j b Bi + (15) 
iso 

where the symbol = indicates isotropization (which is re- 
ally only done for simplicity and should be avoided when it 
becomes important), and dots indicate the presence of further 
terms which here only lead to triple correlation terms. This 
j ■ b term provides an important correction to the usual a ef- 
fect that results from 

(u x b)i = -£ t jkUkUj }P B p + ... = — |w • u Bi + ... (16) 

It is the j ■ b current helicity term whose evolution, under the 
assumption of scale separation, can be described in terms of 
magnetic helicity conservation. Since the magnetic helicity 
is conserved in the large magnetic Reynolds number limit, 
the saturation time of a nonlinear dynamo can be very long. 
This is now described in great detail in several recent reviews 
(Brandenburg 2003, Brandenburg & Subramanian 2004). We 
emphasize however that the j ■ b contribution is quite deci- 
sive in describing correctly the slow saturation phase of any 
nonlinear helical large scale dynamo in a periodic box (Bran- 
denburg 2001, Mininni et al. 2003). 

4. Problem III: a quenching 

4.1. Non-universality of catastrophic quenching 

As shown in Brandenburg (2001), the slow saturation be- 
havior of closed box dynamos can be reasonably well be 
described in a mean field model using catastrophic quench- 
ing, i.e. a(B) = ctoQ(B) and r]%(B) = r)%oQ(B), where 
Q(B) = (1 + Rmlf/B^)- 1 . As discussed in detail by 
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4.2 Dynamical a quenching 



Blackman & Brandenburg (2002), the reason for the agree- 
ment with the simulations is due to the 'force-free degener- 
acy' of the a 2 dynamo in a periodic box, because then B and 
J are then everywhere parallel to each other. This degeneracy 
is lifted for aft dynamos, i.e. in the presence of shear. 

In Fig. [K)] we show the evolution of poloidal and toroidal 
fields, B x and B y , at one point in the simulation of Branden- 
burg et al. (2001), which has shear. Note the systematic phase 
shift and a well-defined amplitude ratio between B x and B y . 
Note also that the dynamo wave is markedly non-harmonic. 
These are clear properties that can be compared with mean- 
field model calculations; see Fig. ^2 A similar type of non- 
harmonic temporal behavior has been found in the first-ever 
nonlinear simulation of a mean field dynamo by Stix (1972); 
see also Fig. ^] An important difference between the two 
models is that for the results shown in Fig.^Jboth a and % 
where quenched, whereas in the model shown in Fig. \l2\m 
was kept constant and a was quenched in a step function-like 
fashion. Since the results from these two very different mod- 
els are similar, the temporal behavior alone cannot really be 
used to discriminate one model over another. 

4.2. Dynamical a quenching 

The a effect formalism provides so far the only workable 
mathematical framework for describing the large scale dy- 
namo action seen in simulations of helically forced turbu- 
lence. The governing equation for the mean magnetic field 
is 



f = Vx(tfxB 



£ - T)(l j) , 



(17) 



where £ = u x b is the electromotive force resulting from 
the u x b nonlinearity in the averaged Ohm's law. MTA can 
be used to derive an expression for £ in terms of the mean 
field. For slow rotation one finds (Radler et al. 2003, see also 
review by Brandenburg & Subramanian 2004) 



£i 



atjBj 



■ VijJj + (7 x B + 5 X J)i + K,./. II,.;.. (18) 

Under isotropic conditions, or as a means of simplification, 
one often just writes 



£ = aB — r/ t J 



(19) 
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Fig. 10. Evolution of B x and B y at x = — tt and z = 0. Note 
that B x has been scaled by a factor —100. [Adapted from 
Brandenburg et al. (2001).] 
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Fig. 11. Evolution of a B B x and 
dimensional mean-field model with a dynamo number T> = 
10, a kinematic growth rate A = 0.015 (which determines 
a) and a microscopic magnetic diffusivity r\ = 5 x 10~ 4 . 
Note that B x has been scaled by a factor 100. (In this case 
the shear S > 0, so we have plotted +B. X , and not — B x as 
we did in Fig.^|where S < 0.) [Adapted from Brandenburg 
et al. (2001).] 




Fig. 12. Evolution of toroidal and poloidal fields, B and A, 
respectively, from the paper by Stix (1972). Note the similar- 
ity in the evolution of his model and that in the direct sim- 
ulation shown in Fig. [K)] with a more spiky behavior of the 
poloidal field compared to the toroidal field. [Adapted from 
Stix (1972).] 

In any case, and this comes as a rather recent realization 
(Field & Blackman 2002, Blackman & Brandenburg 2002, 
Subramanian 2002), there will be an extra magnetic contri- 
bution to the a effect a — > a — qk + om- In the isotropic 
case, aK = — k,T~u) ■ u and ayi = \rj ■ b/p. This leads to 
(Kleeorin & Ruzmaikin 1982; see also Zeldovich et al. 1983; 
and Kleeorin et al. 1995) 

,.2 ( (£ ■ B) t q M ' 



dt 



= -2r) W k 2 



B 2 



(20) 



where we have defined the magnetic Reynolds number as 

Rm = r)ta/r), (21) 

where rj t o = ^t(u 2 ) is the kinematic value of the turbulent 
magnetic diffusivity and B 2 = pqPq(u 2 ) is the equiparti- 



tion field strength. Equation J20i agrees with the correspond- 
ing equation in Kleeorin et al. (1995) if their characteris- 
tic length scale of the turbulent motions at the surface, l s , 
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is identified with 27r/fcf and if their parameter // is identi- 
fied with 87r 2 77j /((tt 2 )/g). The definition of R m may not 
be very practical if r/ t o is not known, but comparison with 
simulations (Blackman & Brandenburg 2002) suggests that 
R m « 0.3i? m , where R m — u Tms /(i~ikf) is a more practical 
definition suitable for simulations of forced turbulence. 

The basic idea is that magnetic helicity conservation must 
be obeyed, but the presence of an a effect leads to magnetic 
helicity of the mean field which has to be balanced by mag- 
netic helicity of the fluctuating field. This magnetic helicity 
of the fluctuating (small scale) field must be of opposite sign 
to that of the mean (large scale) field (e.g. Blackman & Bran- 
denburg 2003). In the following we refer to Eq. J20b as the 
the dynamical a quenching equation (Blackman & Branden- 
burg 2002). Assuming that the large scale magnetic field has 
reached a steady state, and solving this equation for a, yields 
(Kleeorin & Ruzmaikin 1982, Gruzinov & Diamond 1994) 

a K + Vt R m (J-B)/B 2 eq 

a = = (forda/di = 0). (22) 

l + R m (B 2 )/B* q 

Note that for the numerical experiments with an imposed 
large scale field over the scale of the box (Cattaneo & Hughes 
1996), where B is spatially uniform and therefore J = 0, one 
recovers the 'catastrophic' quenching formula, 

a = ^5 (for J = 0), (23) 

l + R m (B )/B c 2 q 

which implies that a becomes quenched when (B ) /B% = 
i?^ 1 10~ 8 for the sun, and for even smaller fields in the 
case of galaxies. 

Obviously, the assumption of a steady state is generally 
not permitted. Especially in the case of closed boxes that are 
so popular for simulation purposes, there is no other way to 
get rid of magnetic helicity than via microscopic resistivity. 
There is now a significant body of literature (see Branden- 
burg 2001 for a comprehensive coverage of the simulation 
results and Brandenburg 2003 for a recent review). To allow 
for faster time scales, and this was realized first by Black- 
man & Field (2000a,b) and Kleeorin et al. (2000), there might 
be hope to achieve this by allowing helicity flux to escape 
through the boundaries of the domain. Two types of results 
will be discussed in the next section. 

5. Open boundaries and shear: the solution? 

In a recent paper, Brandenburg & Sandin (2004) have carried 
out a range of simulations for different values of the mag- 
netic Reynolds number, R m = u Tms /(r]kf), for both open 
and closed boundary conditions using the geometry depicted 
in Fig.^1 In the following we discuss first some results for a 
quenching and an interpretation of the results in terms of the 
current helicity flux. Next we turn to direct simulations of the 
dynamo without an imposed field. 

5.1. Results for a quenching 

In order to measure a, a uniform magnetic field, Bq = const, 
is imposed, and the magnetic field is now written as B = 




x 



B x =B y =B„=0 
A x ,z=A y 2 =A Z =0 

Fig. 13. Differential rotation in our cartesian model, with the 
equator being at the bottom, the surface to the right, the bot- 
tom of the convection zone to the left and mid-latitudes at the 
top. [Adapted from Brandenburg & Sandin (2004).] 

£?o + V x A. Brandenburg & Sandin (2004) determined a 
by measuring the turbulent electromotive force, and hence 
a = (£) ■ Bq/Bq. Similar investigations have been done 
before both for forced turbulence (Cattaneo & Hughes 1996, 
Brandenburg 2001) and for convective turbulence (Branden- 
burg et al. 1990, Ossendrijver et al. 2001). 

As expected, a is negative when the helicity of the forc- 
ing is positive, and a changes sign when the helicity of the 
forcing changes sign. The magnitudes of a are however dif- 
ferent in the two cases: \a\ is larger when the helicity of the 
forcing is negative. In the sun, this corresponds to the sign of 
helicity in the northern hemisphere in the upper parts of the 
convection zone. This is here the relevant case, because the 
differential rotation pattern of the present model also corre- 
sponds to the northern hemisphere. 

There is a striking difference between the cases with 
open and closed boundaries which becomes particularly clear 
when comparing the averaged values of a for different mag- 
netic Reynolds numbers; see Fig. ^] With closed bound- 
aries a tends to zero like i?" 1 , while with open boundaries 
a shows no such decline. There is also a clear difference be- 
tween the cases with and without shear together with open 
boundaries in both cases. In the absence of shear (dotted line 
in Fig. I14> a declines with increasing R m , even though for 
small values of R m it is larger than with shear. The difference 
between open and closed boundaries will now be discussed in 
terms of a current helicity flux through the two open bound- 
aries of the domain. 



5.2. Interpretation in terms of current helicity flux 

It is suggestive to interpret the above results in terms of 
the dynamical a quenching model. However, the quenching 
equation has to be generalized to take the divergence of the 
flux into account. In order to avoid problems with the gauge, 
it is advantageous to work directly with j ■ b instead of a ■ b. 
Using the evolution equation, db/dt = — V x e, for the fluc- 
tuating magnetic field, where e = E — E is the small scale 



5 OPEN BOUNDARIES AND SHEAR: THE SOLUTION? 



5.3 Dynamos with open surfaces and shear 



0.100 



E 

'a 

- 0.010 



0.001 







open 

boundaries 








i a>0 


* 


J>r^ 




N ! 

a<0 . 




-FT 1 ' ^ . 




no shear 






closed 
boundaries 
' 



1 



10 

FL 



100 



Fig. 14. Dependence of |(a)|/u rms on R m for open and 
closed boundaries. The case with open boundaries and nega- 
tive helicity is shown as a dashed line. Note that for R m w 30 
the a effect is about 30 times smaller when the boundaries are 
closed. The dotted line gives the result with open boundaries 
but no shear. The vertical lines indicate the range obtained by 
calculating a using only the first and second half of the time 
interval. [Adapted from Brandenburg & Sandin (2004).] 



electric field and E = tj^qJ — S the mean electric field, one 
can derive the equation 



^-j b = -2e—c-V ■ Ff?, 
dt ° 

where 



C 



2e x j + (V x e) x b/fi , 



(24) 



(25) 



is the current helicity flux from the small scale (SS) field, 
and c = V x j the curl of the small scale current density, 
j = J — J. In the isotropic case, e • c w kf e ■ b, where k{ is 
the typical wavenumber of the fluctuations, here assumed to 
be the forcing wavenumber. 

Making use of the adiabatic approximation one arrives at 
the algebraic steady state quenching formula (da/dt = 0) 



a K + R m (rjtVoJ ■ B - \k t 2 V • Mo^c) /&> 



1 + R m B 2 /B2 



(26) 



In the absence of a mean current, e.g. if the mean field is 
defined as an average over the whole box, then B = Bq = 
const, and J = 0, so Eq. J26l > reduces to 



a K -\ki 2 R m V ■ n Tf/B? 



cq 



1 + R ni B /B% 



(27) 



This expression applies to the present case, because we con- 
sider only the statistically steady state and we also define the 
mean field as a volume average. 

In the simulations, the current helicity flux is found to be 
independent of the magnetic Reynolds number. This explains 
why the a effect no longer shows the catastrophic i?" 1 de- 
pendence (see Fig. I14> . In principle it is even conceivable 
that with ax = a current helicity flux can be generated, for 
example by shear, and that this flux divergence could drive a 
dynamo, as was suggested by Vishniac & Cho (2001). It is 
clear, however, that for finite values of R m this would be a 




Fig. 15. Visualization of the toroidal magnetic field during 
three different times during the growth and saturation for the 
run without kinetic helicity. 



non-kinematic effect requiring the presence of an already fi- 

l li 

nite field (at least of the order of B cq / R^ ). This is because 
of the 1 + i? m Bg/ B^ q term in the denominator of Eq. i27l . 
At the moment we cannot say whether this is perhaps the ef- 
fect leading to the nonhelically forced turbulent dynamo dis- 
cussed in Sect. 15.31 or whether it is perhaps the d x J or 
shear-current effect that was also mentioned in that section. 

5.3. Dynamos with open surfaces and shear 

The presence of an outer surface is in many respects simi- 
lar to the presence of an equator. In both cases one expects 
magnetic and current helicity fluxes via the divergence term. 
A particularly instructive system is helical turbulence in an 
infinitely extended horizontal slab with stress-free boundary 
conditions and a vertical field condition, i.e. 



l y,z — u z — B x — B y — 0. 



(28) 



Without shear, such simulations have been performed by 
Brandenburg & Dobler (2001) who found that a mean mag- 
netic field is generated, similar to the case with periodic 
boundary conditions, but that the energy of the mean mag- 

2 

netic field, (B ), decreases with magnetic Reynolds num- 
ber. Nevertheless, the energy of the total magnetic field, 
(B 2 ), does not decrease with increasing magnetic Reynolds 
number. Although they found that (B ) decreases only like 

— 1/2 

R m , new simulations confirm that a proper scaling regime 
has not yet been reached and that the current data may well 
be compatible with an R m l dependence. 

Clearly, an asymptotic decrease of the mean magnetic 
field must mean that the small scale dynamo does not work 
with such boundary conditions. Thus, the anticipated advan- 
tages of open boundary conditions are not borne out by this 
type of simulations. In the presence of shear the results are 
very different. We now use the same setup as in Sect. 15.11 
The size of the computational domain is \ix x 2-7T x i-7r and 
the numerical resolution is 128 x 512 x 128 meshpoints. The 
magnetic Reynolds number based on the forcing wavenumber 
and the turbulent flow is around 80 and shear flow velocity 
exceeds the rms turbulent velocity by a factor of about 5. We 
have carried out experiments with no helicity in the forcing 
(labeled as a = 0), as well as positive and negative helicity 
in the forcing (labeled a < and a > 0, respectively); see 
Fig. for a visualization of the run without kinetic helic- 
ity. We emphasize that no explicit a effect has been invoked. 
The labeling just reflects the fact that, in isotropic turbulence, 
negative kinetic helicity (as in the northern hemisphere of a 
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Fig. 16. Saturation behavior of the ratio q — (B")/(B 2 ) 
for runs with different kinetic helicity of the flow. Solid line: 
zero helicity, dotted line: positive helicity (opposite to the 
sun) dashed line: negative helicity (as in the sun). The line 
at denoted by "a = 0, closed" refers to a case where the nor- 
mal field condition on the equator and the surface has been 
replaced by a perfect conductor condition. 




Fig. 17. Gray/color scale representation of the azimuthally 
and time averaged mean azimuthal field B together with vec- 
tors in the meridional plane. In this run the turbulence is non- 
helically forced (a = 0). 



star or the upper disc plane in galaxies) leads to a positive a 
effect, and vice versa. 

We characterize the relative strength of the mean field by 
the ratio q = (B )/(B 2 }, where overbars denote an average 
in the toroidal (y) direction; see Fig.^] A strong mean field 
is generated in all cases, unless a perfect conductor boundary 




Fig. 18. Same as Fig.^J but for a > 0, so the turbulence is 
forced with negative helicity. Note the clockwise sense of the 
poloidal field. Together with the positive toroidal field (indi- 
cated by light shades) this corresponds to positive magnetic 
helicity, consistent with a > 0. 




Fig. 19. Same as Fig. El but for a < 0, so the turbulence 
is forced with positive helicity. Note the anti-clockwise sense 
of the poloidal field, corresponding to negative magnetic he- 
licity, consistent with a > 0. 



condition ("closed") is adopted on the outer surface and on 
the equator. The mean field appears to be statistically station- 
ary in all cases, i.e. there is no indication of migration in the 
meridional plane. A time average of the mean field is shown 
inFig.EI 

There are two surprising results emerging from this work. 
First, in the presence of shear rather strong mean fields can 
be generated, where up to 70% of the energy can be in the 
mean field; see Fig.^] Second, even without any kinetic he- 
licity in the flow there is strong large scale field generation. 
However, this cannot be an afl dynamo in the usual sense. 



6 CONCLUSIONS 




Fig. 20. Vectors of <fii yy together with contours of U y which 
also coincide with the streamlines of the mean vorticity field 
W . Note the close agreement between <f> iyy and W. The ori- 
entation of the vectors indicates that negative current helicity 
leaves the system at the outer surface (x = 0). 



One possibility is the 8 x J effect, which emerged originally 
in the presence of the Coriolis force; see Radler (1969) and 
Krause & Radler (1980). In the present case with no Coriolis 
force, however, a 6 x J effect is possible even in the presence 
of shear alone, because the vorticity associated with the shear 
contributes directly to£ocVF = Vx£/ (Rogachevskii & 
Kleeorin 2003, 2004). 

A rather likely candidate for the current helicity flux is the 
so-called Vishniac & Cho (2001) flux. A systematic deriva- 
tion using MTA (Subramanian & Brandenburg 2004) gave 

GQ 

Fc^fakBjBk, (29) 
where tfiijk is a new turbulent transport tensor with 



<j>ijk = -4TU>kVjUi 



(30) 



Obviously, only components of that are symmetric in j 
and k enter the flux ^F(f . Therefore we consider in the fol- 
lowing <f>y k = 7}(4>ijk + 4>ikj)- In the present case where the 
toroidal field (y direction) is strongest, the <fii yy components 
(where i = x and z) are expected to be most important. The 
time average of this component is shown in Fig.|20| Here we 
have only considered the nonhelical case, but the helical cases 
are indistinguishable within error margins. All the 6 indepen- 
dent components of cj^ij and ^ - fe are given by 



Note that the trace of both tensors is small. Thus there is 
no hope that this tensor could possibly be isotropic. This is 
because isotropization, i.e. ^8jk<f>ijk> would then give some- 
thing close to zero. Nevertheless, one may hope that some 
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useful parameterization of (pfjf., that can be used in mean field 
calculations, will soon be available. 

In conclusion, there is evidence that the strong dynamo 
action seen in the simulations is only possible due to the com- 
bined presence of open boundaries and shear. 



6. Conclusions 

In the present work we have discussed three aspects that are 
strongly connected with the applicability of mean field dy- 
namo theory. The first concerns turbulent diffusion, i.e. the 
ability of the flow to mix a large scale magnetic field such 
that its energy can be converted into energy at smaller scales. 
In order that microscopic diffusion can eventually be ther- 
malized, unimpeded by the Lorentz force, the field has to 
be weak enough at the scale where diffusion takes place. 
Whether this is indeed the case can in principle be seen from 
the simulations. Since in the simulations the magnetic and 
fluid Reynolds numbers are much smaller than in reality, it 
would be important to show that the spectra begin to con- 
verge for large magnetic and fluid Reynolds numbers. At the 
moment this is simply not yet the case, even at a resolution 
of 1024 3 meshpoints. However, by a combination between 
direct and large eddy simulations some preliminary insight 
can be gained. The suggestions are that magnetic and kinetic 
energy spectra approach spectral equipartition in the deep in- 
ertial range, and that the slope would be comparable to the 
Kolmogorov k~ 5/>3 slope. 

The other aspect discussed in this paper concerns the cal- 
culation of turbulent transport coefficients such as the a ef- 
fect. It has been particularly disturbing that much of the the- 
ory in this field was based on the first order smoothing ap- 
proximation which clearly should break down under physi- 
cally interesting conditions. The reason for its apparent suc- 
cess is now becoming clear and are connected with the close 
similarity between the first order smoothing approximation 
and the minimal tau approximation where the triple correla- 
tions are no longer omitted. The reason for the close simi- 
larity between the two approaches is mainly connected with 
the fact that all the interesting physics enters via the quadratic 
correlations which can be worked out correctly using linear 
theory. The significance of the triple correlations is merely 
that they determine the length of r. In the framework of the 
first order smoothing approximation the value of r had to 
be obtained by dimensional arguments and was therefore not 
well determined. 

Next we have briefly addressed the long standing issue 
of a quenching and have pointed out the connection be- 
tween dynamical quenching and catastrophic quenching. The 
possibility of catastrophic a quenching is closely linked to 
magnetic helicity conservation. The degree of quenching can 
therefore be alleviated by allowing for magnetic or better cur- 
rent helicity fluxes. Such fluxes do not just come on their own, 
so simply allowing for the boundaries to be open is not suf- 
ficient. One needs to drive a flux through the entire domain. 
This seems to be accomplished by the Vishniac-Cho flux. As 
already argued by Vishniac & Cho (2001), and confirmed by 
Arlt & Brandenburg (2001), this flux can be driven in the 
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presence of shear. Indeed, there are strong indications that 
the helicity flux can alleviate catastrophic quenching at least 
by a factor of about 30, and that it allows for dynamo ac- 
tion generating strong magnetic fields, even in the absence of 
helicity in the forcing. The details of this process are not en- 
tirely clear, but possible candidates include the shear-current 
effect of Rogachevskii & Kleeorin (2003, 2004), and perhaps 
the helicity flux itself (Vishniac & Cho 2001) which is a non- 
linear effect. 

In the near future it should be possible to investigate the 
emergence of current helicity fluxes from a dynamo simu- 
lation in more detail. This would be particularly interesting 
in view of the many observations of coronal mass ejections 
that are known to be associated with significant losses of 
magnetic helicity and hence also of current helicity (Berger 
& Ruzmaikin 2000, DeVore 2000, Chae 2000, Low 2001, 
Demoulin et al. 2002, Gibson et al. 2002). In order to be able 
to model coronal mass ejections it should be particularly im- 
portant to relax the restrictions imposed by the vertical field 
conditions employed in the simulations of Brandenburg & 
Sandin (2004). One possibility would be to include a sim- 
plified version of a corona with enhanced temperature and 
hence decreased density, giving rise to a low-beta plasma ex- 
terior where nearly force-free fields can develop (Gudiksen 
& Nordlund 2003). Such a setup might allow detailed com- 
parison with observations. 
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